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ON ■ Abstract 

Q\ ■ 

We utilize asymptotic Pade-approximant methods to estimate the three-loop order MS-scheme coefficients 

i-Q ' within the inclusive b — > uvif decay rate for four and five active quark flavours. The estimates we obtain 

^Hi for the three renormalization-group-accessible coefficients within the three-loop contribution are all found to 

be within 5.1% of their true values, using a least-squares procedure in conjunction with an asymptotic Pade- 

Q_) \ approximant estimate of the three-loop term over the entire [i > mi, domain. Given the input values a s (Mz) — 

,_J^ ■ 0.118±0.004 and mb{mb) = 4.17±0.05 GeV, the three-loop expression for the purely-perturbative contribution 

to the b — > ui>e£~ decay rate is minimally sensitive to renormalization scale at [s, — 1.775 GeV, at which scale 

the three-loop contribution is estimated to be only 1.4% of the leading tree-order contribution. We estimate the 

full perturbative decay rate to be 1927r 3 r(6 -» uvil~)/ (G 2 F \V ub \ 2 ) = 2070 GeV 5 ± 16%, inclusive of theoretical 

uncertainties from series truncation, the input parameters, and the estimation procedure. 



One of the outstanding issues in B physics is the accurate determination of the CKM matrix element V u b, whose 
magnitude corresponds to one side of the unitarity triangle. Experimental determination of the angles within this 
triangle is anticipated soon from i3-factories, i.e. BELLE and BaBar detectors, and from experiments by the CDF 
and CLEO III collaborations. A more precise determination of \V u b\ can test the consistency of of the CKM picture 
of quark mixing within the standard model. The value of \V u b\ can be extracted from the inclusive semileptonic 
width T(B — ► X u v^) once a more accurate determination of this decay rate becomes available. In order to 
reduce the theoretical uncertainty in this rate, it is useful to determine its purely perturbative QCD corrections as 
accurately as possible. 

In perturbative QCD, the inclusive semileptonic b — > u decay rate T{b — > uvg(.~) may be expressed as pi 



— r 0, m b {p),x(fi)) = m\{n) ( 1 + [oo - ai log(w)] x{p) + [b - h log{w) + b 2 log 2 (w)] x 2 (n) 



+ [c - Cl log{w) + c 2 log 2 {w)-c 3 log 3 (to)] x 3 (p) + . . . j 



(1) 
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where 

_ ocjjj) , , ,, _ ml(fi) _ G\ \V ub \ 2 

X W - — — ' w = w{n,m h {n)) = — — , K = i927r3 . (2) 

If mj is identified with a scale-independent pole mass, the known first and second order terms within (pi) appear to 
be poorly convergent because of the proximity of a Borel-plane singularity 0, 0| ■ In ref . M , this problem is averted 
by employing MS renormalization in which mi, is identified with a scale-dependent (running) 6-quark mass for four 
or five active flavours. The one- and two-loop order MS coefficients within (fil) are then given by fm 



all Uf : 


a = 4.25360 , 


oi = 5 , 






71/ = 5 : 


6o = 26.7848 


, 6i = 36.9902 


b 2 ~- 


= 17.2917 


n/ = 4 : 


6o = 25.7547 


6i = 38.3935 


b 2 -- 


= 17.7083 



(3) 

We utilize asymptotic Pade-approximant methods in the present note to extract an estimate of the three- 
loop order contributions to (fy). Such estimates are then tested against explicit (MS) renormalization-group (RG) 
determinations of c\, c%, and C3 (the coefficient cq is RG-inaccessible to order a; 3 ). We also find, as in M, that the 
optimal value of the scale parameter \i (determined by scale-insensitivity) for evaluating (nl) is substantially below 
rrib(mb) [defined to be the solution to fi = rrib(p)]. We then present an estimate of the magnitude and theoretical 
uncertainties of the inclusive semileptonic b — > u decay rate ([l]). 

To begin, consider a truncated perturbative series 



JV+1 



S N+1 = 1 + > R n x n , S = lim S N (4) 

n=l 

in which {i?i, R2, ■ ■ ■ , Rn+i} are known. We seek to utilize Pade approximant methods to estimate the first 
unknown term Rjy + 2X N+2 within the series S. For example, the set of coefficients {R\, R2, ■ ■ ■ , Rn+i} is sufficient 
to generate an [N\l] approximant to the series S, 

l + a 1 x + ... + a N x N 

S[Nm = TTm ' () 

such that the first TV + 2 terms of the Maclaurin expansion for S[nw correspond to the truncated series Sjv+i- The 
next term of this Maclaurin expansion represents a Pade-approximant estimate of Rn+2% N+2 ■ Specifically, we find 
that 

p [JV|i] _ ■gjy+i ,„s 

h n+2 - S — ■ ( b > 



Rn 

To improve such an estimate, we note that the relative error of an [7VJM 
in a perturbative field-theoretical series is anticipated to be || [|, [|, f| ^ 



-approximant estimate [as in (ph] to terms 



n-Tjj-M+i _ -fw+M+i __ MIA 

R N+M +i (N + M(l + a) + b) 1 



c[N\M] __ f%+M+l ~ -tW+M+l _ ivL-.n 

°N+M+1 ~ p. -~/»r , w/, , s . rxAf ' V'J 



In (0), i?]y ! M+1 is the estimate for i?7v+A/+i obtained from a Maclaurin expansion of the [iV|M] approximant 
generated by the first N + M terms of a series S. Consequently, 5^+m+i is the relative error associated with such 
a Pade estimate, and the constants A, a, and b serve to parameterize this relative error. We now use (0) to obtain 
an improved estimate of the three-loop coefficient R3 in a series, given knowledge only of the one- and two-loop 
coefficients R\ and R2. If M — 1 and if we assume (as in prior applications |7], ||) that a + b = 0, we see from 
substitution of (g) into (0) that 

4 0|11 = ^^ = -^ , (8) 

s li\i ] = l%/Ri-R 3 = _ A/2 (9) 

Ra 



Substituting (g) into (g), we see that 

R 3 



2R\ 



R l {R\ + R 2 ) 



(10) 



For the perturbative series on the right-hand side of (fy), both the known coefficients i?i i2 and the unknown 
coefficient i? 3 are explicitly of the form [see (ffl)]: 



Ri{w) = a a -ailog(iu) , 

R 2 (w) = 6 -&ilog(iu) +6 2 log 2 (w) , 

•R3M = c - ci log(» + c 2 log 2 (u;) - c 3 log 3 (w) 



(11) 
(12) 

(13) 



To obtain the unknown coefficients Co- 3 within the three- loop order term R^(w), we compare moments of Rsiw) 
over the ultraviolet region < w < 1 [i.e., fi> mj(m(,)], as obtained from substitution of ( p"T| ) and (|l^) into (|J0|), 
to the same moments obtained explicitly from (|13l): 



N k = (k + 2) I dw w" 





fe+i 



2^i(w) 



_i?iH(i? 2 (u)) + i? 2 (w)). 

= (fc + 2) / dw w k+1 [c - ci log(w) + c 2 log 2 (w) - c 3 log 3 (w)] 



(14) 

Using the nj = 5 values given in (0) for the known coefficients {ao-i, 60-2}, we obtain the following set of equations: 
7V_i = 1118.64= c + ci + 2c 2 + 6c 3 , (15) 

(16) 

(17) 

(18) 



113 
N = 457.181 = c + -a + -C2 + -cs 

JVi - 337.233 = c + - Cl + -c 2 + -c 3 

1 1 3 

N 2 = 291.645 = c + ~c x + -c 2 + — c 3 
4 o oz 



The solution of these equations are our estimates of the three-loop coefficients in dl]) for the case of five active 
flavours: 



i f = 5(est) : c = 200.5 , a = 251.4 , c 2 = 190.5 , c 3 = 47.61 



(19) 



True values for the coefficients ci_ 3 (but not Co) can be extracted from the invariance of the physical transition 
rate (||) under changes of scale \i JO , 



= fi—T (n, m b (fi), x(n)) 
d/j, 



d , 2 , s d 

Mo (7 a; + 7ix +722: + . . .) m b - — 

O/i orrib 



(Pox 2 + ftx 3 + . . .) 



d_ 
dx 



r ,(20) 



(21) 



an equation whose order-x 3 contributions vanish provided 

c\ = 26 /3o + a ft + 70 (5&o - 2&i) + 71 (5a - 2oi) + 572 , 
C2 = - [26i/3 + 01 (ft + 57i) + 7o (56i - 46 2 )] , 

c 3 = ~ (2ft + 5 7 o) • 

Using the n f = 5 values in (§) and the values ft = 23/12, ft = 29/12, 7o = 1, 7l = 253/72, and 72 = 7.41985 for 
five active flavours, we find from (Pll) that 



n/ = 5 : ci = 249.592 , c 2 = 178.755 , c 3 = 50.9144 



(22) 



The estimates (19[) for ci_3 are all within 7% of these correct values.R 

Since we do not know the true value of Cq, it is important to test whether the value obtained in (119) is stable 



when the estimates of Ci_3 derived from (|15|— 118|) are replaced with the true values given in (E2). If we utilize (22) 
explicitly within (|J"5|-|l8|) , we obtain four separate determinations of Cq, which are respectively given by c = 206.0, 
204.8, 203.0 and 202.1. These results are all within 1% of each other, and 1-3% above the estimate obtained in 



(19). A further consistency check on the cq estimate is to construct a least-squares fit of equation (pXl|) to the form 



of eq. (O) over the entire region (i > mb(rrib): 



X (co,ci,c 2 ,c 3 ) = / dw 

(i 



1 2 

(c - ci log(u;) + c 2 log (w) - c 3 log {w)j 



[Rt(w) (Rl{w) + R 2 {w)) 
The minimization requirements d\ 2 jdci = lead to the following estimates for c,: 



(23) 



if = 5(est) : c = 198.4 , a = 260.6 , c 2 = 183.9 , c 3 = 48.64 , xL™ = L505 ( 24 ) 



Note that ci_3 are all within 4.5% of their true values (123), and that the small value for Xmin implies a precise 



fit of ( p. OP to the form of (|13| ) over the entire interval < w < 1. The estimates of cq obtained in (]24j) and (19) 
are remarkably close. If true values ( |22| ) for ci_3 are substituted directly into x 2 , minimization with respect to the 
only remaining variable cq leads to the estimate 

n f = 5{est.) : c = 206.0 , (25) 

identical to that obtained by a similar substitution into (|15|). We shall focus on this particular estimate when 
estimating the rate (0) , despite the fact that it is 3-4% above the estimates in (|l^) and (|2J) , because this estimate 
incorporates the RG-determinations of values for ci_3. 

The procedures delineated above can also be applied to the case of four active flavours. Utilization of the nf = 4 



values for {ao-i, ^0-2}, as tabulated in (H), in the first integrand of (14) leads to the following numerical values 
for the moments Nk'. 

AT_ X = 1165.14 , 7V = 457.864, N x = 329.529 , N 2 = 280.696 . (26) 

If these values are incorporated into (|l5|-|l8|) , the solution to the resulting set of equations is 

n f = 4(est.) : c Q = 182.7 , a = 271.9 , c 2 = 201.4 , c 3 = 51.26 (27) 

True values for ci_3 are RG-accessible via eqs. (pi|), with nf = 4 values for {ao-i, ^0-2} as listed in (g), and with 
(3q = 25/12, j3\ = 77/24, 70 = 1, 71 = 263/72, and 7 2 = 9.94702 for four active flavours. The correct values are 

n f = 4 : a = 263.84 , c 2 = 194.23 , c 3 = 54.109 (28) 



which are within 5.3% of their estimated values in (|27|). As before, the incorporation of the correct values (|28|) into 
the moment equations (jl5]-|lg|) with numerical values ( pq ) for rif = 4 moments leads to four separate estimates for 
Co {188.2, 188.3, 186.4, 185.4} that are remarkably consistent, and all less than 3% above the estimate in (p7[). 
Moreover, minimization of the \ 2 function ( p3[ ) using Jo) and ( |12| ) with nf = 4 values for {ao-i, &o- 2 } yields 
estimates for ci_3 all within 5.1% of their true values (|28|); i.e., d\ 2 /dci — when 

n f = 4(est.) : c = 181.5 , a = 277.3 , c 2 = 197.6 , c 3 = 51.86 , xLn = 0.5079 (29) 



The above estimate for Co is within 1% of that obtained in (^7|). If one incorporates into % 2 the true values ( J28[ ) 
for Ci_3, one finds that d\ 2 /dco — when 

nf = 4(est.) : c Q = 188.2 , (30) 



1 The moment methodology delineated above has also been employed with similar success to the RG-accessible order-x 3 terms of 
the gluonic scalar-current correlation function [hcj and the CP-odd Higgs boson decay rate [hl|. 



analogous to the estimate (Eq) for five active flavours. 

To examine the consistency of the Pade-approximant methodology leading to (p9), one can obtain an estimate 
of the rij — 4 value for Cq directly from ( p5| ) via comparison of four- and five-active-flavour expressions for the decay 
rate (El). Suppose for n/ = 5 [indicated henceforth by the superscript ' 5 '] that we evaluate (fil) at \i — 4.17 GeV in 
order to coincide with (the central value of) an empirical estimate of m b (nib) |Q, and we choose a;^ 5 ^(4.17GeV) = 
0.0715492, a value devolving (numerically) via the four-loop /3-function from a s (Mz) — 0.118 (the corresponding 
input value in fL2| ), If we substitute into (p the n/ = 5 estimate ( |25|) for co and known n/ = 5 values (0) for ao 
and bo, and if we note that log(u>) = when /i = m b (/x), the defining relationship for m b (nib), we find that 



rO) 

-— = (4.17 GeV) 5 
K 



1 + 0.304342 + 0.137119 + 0.075454] = 1912.68 GeV 5 



(31) 



The same rate can be evaluated for the case of four active flavours with cq left arbitrary by substituting into 
P) the rif = 4 RG values (£8) for ci_3 and the rif = 4 values (0) for {ao_i, &0-2}, and by utilizing following 
threshold-matching conditions Il3[ to obtain values for x(4.17GeV) and m(,(4.17GeV) appropriate for four active 
flavours, as indicated by the superscript ^' : 



,(4) 



,(4) 



(4.17 GeV) = 



(4.17 GeV) = 



= ~(5) 



(4.17 GeV) 



r t2 r -1 3 

1 + 0.1528 x (5) (4.17 GeV) +0.633 x {5) (4.17 GeV) 



0.0716218 , 
m ( b ] (4.17 GeV) 
4.17739 GeV 



1 + 0.2060 



,(5) 



(4.17 GeV) 



1.9464 



,(5) 



(4.17 GeV) 



(32) 



(33) 



We note that log(iy) = 2 log 
expression for the rate: 



,(4), 



(h)//j, is no longer zero, and we obtain from ([l]) the following four-active-flavours 



r(4) 

— = (1824. 



0.46737c ) GeV 5 



(34) 



Equating (31) to (B4h, we find that cq in the latter expression is 188.03, a value in startlingly close agreement 



,(5), 



with the independent estimate (|30|). Thus, at the five-flavour threshold ml (rrib) = 4.17GeV, equilibration of the 
below-threshold and above-threshold expressions for the inclusive semileptonic b — > u rate is seen to corroborate 
the consistency of the independent cq estimates ( |30| ) and (|2q ) for four and five active flavours. 

As noted in jl[], the series (|]) will exhibit some u, dependence as an artifact of its truncation; i.e., the three 
loop expression for T is subject to 0(x 4 ) violations of (pp|). In Figure 1, we display this residual fi dependence 
of the rate T[u]/K both below and above the five-flavour threshold, which is assumed (as above) to occur at 



,( 5 ) 



(m b ) = 4.17 GeV. For fi greater than this threshold, we obtain T\p\/K utilizing known nt = 5 values 



for 



{ao-i, &o-2j £1-3} hi conjunction with the estimate (|25| ) for c . The above-threshold values for m b (jj,) and x^'(/j) 
are evolved numerically via rif — 5 four-loop RG-functions from the specific values m h (4. 17 GeV) = 4. 17 GeV 



'% v 
and x^(Mz) = 0.118/7T. For u, less than 4.17GeV, we utilize known nf = 4 values (| 28) for {ao-i, 60-2, C1-3} 

in conjunction with the estimate ( p0| ) for cq. The below-threshold values for x^(u) and m b (p) are evolved via 



,(4) 



rif = 4 four-loop RG-functions from the specific values for x 

(33). We note that the maximum of the rate exhibited in Fig. 1 occurs at [i 

that 



(4.17 GeV) and m[ 4) (4.17 GeV) obtained in 



% and 



T [1.775 GeV] 
K 



[5.1213 GeV] 5 [1 - 0.6455 + 0.2477 - 0.0143] = 2071 GeV 5 



1.775 GeV (near the r mass), and 



(35) 



This value corresponds to the minimal-sensitivity Q prediction of the decay rate. 

In order to get a handle on the theoretical uncertainty of this estimate, we will assume that there exists a bi- 
directional series truncation error equal to the three- loop contribution to (|3^). An additional source of theoretical 



uncertainty is the error associated with the estimate (BOh for cq. This uncertainty should be comparable to the 
explicit error in estimating c\, C2, and C3 by the same least squares fitting procedure in (p9[). Comparing the true 



values (28) to these estimates, we have already noted that such errors are 5% or less, a value we shall assume to 
characterize error in the estimate (pQ). Such an error estimate is upheld by the range 198 < cq < 206 obtained 
in the different methodological approaches to determining cq [see (pl]j2q)]. As noted above, the estimate (35) 



is based upon (correlated |12]|) input values a s (Mz) = 0.118 and m b (nib) = 4.17 GeV. The errors associated 



with these estimates are ±0.004 and ±0.05 GeV, respectively [12(] . Finally, we quote the ref. |15j estimate of the 
nonperturbative (NP) contribution to T / K: 

*£\ =m f f -9*2 + Ai \ ^ _o.5GeV 2 <A!<0 , A 2 = 0.12 GeV 2 (36) 

K J NP V 2 / 

Taking into account all of these uncertainties, we estimate that 

im^nb^uvd-) = 2Q71 GeV 5 ± ^ GeV 5 ± 35 QeV 5 (+144) GeV 5 (+120) GeV 5 (-39) Gey 5 (37) 

Gp \Vub\ 

where the listed set of uncertainties respectively correspond to uncertainties associated with series truncation, the 



Cq estimate, the value for a s (Mz), the value for rn^^nib), and the nonperturbative contribution (36) 



The truncation error listed in ( p7|) may well be excessive. The series in (|35|) appears to be an alternating series 
whose final term is quite small. If this alternating character is retained for higher orders of x, then the series within 



(35) is bounded from above by the sum of its first three terms and bounded from below by the sum of all four 



known terms, in which case the truncation error should be unidirectional: 

2071GeV 5 <( — j <2122GeV 5 (38) 

\ -**■ / true 

The above suggests that truncation error at this point of minimum scale-sensitivity may lead to a spread as small 
as ±1.3%. However, the smallness of the order- a; 3 contribution suggests that the series ( |35| ) may be asymptotic, 
consistent with the anticipated behaviour (Rn ~ N\B~ N N C |16[) underlying eq. <Q). If such is the case, the 
magnitude of the order-x 4 -and-higher terms may well increase with order. Since the value for a quantity represented 
by an asymptotic series is obtained by a sum of only its decreasing terms (with an error comparable to the smallest 
such term), an optimal determination of (T/K) true may well be obtained by choosing fj, such that the order-a: 3 term 
vanishes. This value of fi = 1.835 GeV corresponds to the point of Figure 2 at which the two-loop and three-loop 
expressions for T/K cross: 

r(1.835GeV) =2069GeV5 
K 

The estimates (pq) and (139) are 8% larger than the rate (pfi) at the five-flavour threshold. The former values 
represent a superior estimate for the true rate, since the rate at this point is locally insensitive to the renormalization 
scale parameter u,. The discrepancy between (pW and (Bq) can be understood by noting that ( J3l| ) appears to be a 
positive term series. If convergent, a truncated positive term series is necessarily an underestimate of the series sum. 
Indeed, the rate displayed in Fig. 1 falls off even more (as noted in H ) for values of /j, above the five-flavour threshold 
m b (nib). This progressive fall off can be understood by noting that the entire set of coefficients {ao-i, 60-2, C0-3} 
is positive, and that when u, increases above the five-flavour threshold, \og(w) becomes increasingly negative. 
Consequently, we see from (II]) that the coefficients of x, x 2 , and a; 3 become increasingly positive, more than 
offsetting the decrease in x itself with increasing u,. This serves to increase the magnitude of the successive terms 
within the positive term series, which necessarily implies an increase in the truncation error, as well as a progressive 
underestimation of the full rate, as a consequence of truncation to three-loop order. 

One way to obviate these difficulties and to recover (approximate) scale invariance of the full series for the rate 
r is to estimate this sum via the use of Pade approximants. For example, such Pade-summation of the known 
series within the Bjorken sum rule has been shown j|] to reduce substantially the scale dependence of the sum- 
rule-extracted value for a s . This reduction in scale dependence via the substitution of (diagonal or near-diagonal) 



Pade approximants in place of truncated field-theoretical perturbative series has been more generally established 
in reference H. 

For the case of the inclusive semileptonic b — > u decay rate, we display in Figure 3 the three-loop truncated 
rate together with its corresponding [1|2] and [2|1] approximants [i.e., the [1|2] and [2|1] Pade-summations of the 
truncated perturbative series S3, as defined by (El)]. The approximants are seen to exhibit substantially less in- 
dependence than the three loop expression for the rate, and are seen to be reasonably close to the minimal-sensitivity 
value obtained in (p5|). The [2|1] approximant to the rate remains between 2001 GeV 5 and 2029 GeV 5 over the 
range 4.17 GeV < /1 < 9 GeV, and the [1|2] approximant remains between 1956 GeV 5 and 2053 GeV 5 over the same 
range. By contrast, the truncated series itself is seen from Figure 3 to fall from 1913 GeV 5 at fi = 4.17 GeV to 
1733 GeV 5 when (j, = 9 GeV. 

At present, the experimental uncertainties in direct measurements of the inclusive semileptonic b —* u decay 
rate are substantially larger than the theoretical uncertainties listed in (37); e.g., recent ALEPH |17j and L3 
flqj Collaboration branching-ratio measurements are (1.73 ± 0.78) • 10 -3 and (3.3 ± 2.0) • 10~ 3 , respectively, with 
systematic and statistical errors combined via quadrature. Somewhat less direct measurements based upon endpoint 
region measurements of the T(4S) decay spectrum are characterized by additional model dependence, although 
they are utilised as a source of present experimental bounds [[19) on V u b : |K,f,/V c b| = 0.08 ± 0.02 and \V u b\ = 
0.0037 ± 0.0006. However, as experimental uncertainties in the semileptonic b — > u rate are reduced, control over 
the uncertainty within theoretical expressions for the rate becomes increasingly important .n The estimate we 
present above for the unknown three-loop coefficient cq is corroborated both by concomitant success in predicting 
the RG-accessible three- loop terms within the theoretical rate (1), as well as by consistency between the direct 
four-active-flavour (asymptotic-Pade) determination of Cq (30) and its independent determination via (25) and (34) 
[i.e. via flavour-threshold matching conditions]. Given the total b-decay rate T tota i — (4.249 ± 0.039) • 10~ 13 GeV 
KM , our result (37) leads to the following relationship between \V u b\ and the inclusive semileptonic branching ratio 
(BR): 

RR-IT/ v( G 2 f \ [(2071 pert . -48 nonpert .)GeV 5 ±l7%] 
1 ubl \192i^J 4.25 • 10- 13 GeV 

\V ub \ = 0.0959 • (BR) 1/2 ± 9%. (40) 
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Figure 1: The /i-dependence of the three-loop order decay rate T(b — » uPe£~)/K obtained from input values 
a s (Mz) = 0.118, mb(mb) = 4.17 GcV. The rate is calculated using four active flavours below and five active 
flavours above the five-flavour threshold, subject to threshold matching conditions discussed in the text. All scales 
are in GeV units. 
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Figure 2: The two-loop (2L) and three-loop (3L) expressions for T(b — > we,l~)/K obtained from input values 
a s (Mz) = 0.118, mb(mb) — 4.17 GcV in the region of minimal ^-sensitivity of the three-loop rate. Both curves are 
obtained using four active flavours. All scales are in GeV units. 
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Fig.3 
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Figure 3: Comparison of three-loop (3L) and Pade-summation expressions for the rate T(b — > uD(£~)/K obtained 
from input values a s (Mz) — 0.118, mb(rrib) = 4.17 GeV when fj, is above the five-flavour threshold. The summations 
are [1|2] and [2|1] approximants whose Maclaurin expansions reproduce the 3L expression. All scales are in GeV 
units. 
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